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Abstract 

Qs^ • The problem of two nonrelativistic chromoelectric and chromomag- 

0^ . netic charges in a Higgs vacuum is considered in classical field theory. 

An approximation of interaction potential is constructed on the basis 
of a numerical solution to the equations of gluostatics. The concept 
of a non-Abelian Abrikosov vortex is discussed. It is shown that the 
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results of Ginzburg-Landau theory for the tension of a string between 
magnetic charges can be directly extended to the non-Abelian case. 
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S Introduction 

OO I At present, the mechanism of dual superconductivity is one of the most 
appealing concepts that are invoked in attempts at explaining quark confine- 
■ ment [|I|. It is assumed that the physical vacuum hinders the penetration of 
gluon fields to large distances from the sources of charge and compresses lines 
Ph! of force in the same manner as tube does. This pattern is well known in the 
QhI Abelian theory of superconductivity for magnetic charges. The mean-field 
^ \ approximation as exemplified by Ginzburg-Landau theory is quite sufficient 
for obtaining a quantitative description of arising physical situation-and in 
particular, for evaluating the tension coefficient 0. Since the microscopic 
theory of confinement has yet to be developed, the above furnishes sufficient 
5^ I motivation for performing a meanfield analysis of the problem of non-Abelian 
charges in a vacuum that possesses the properties of a gluon superconductor. 
In this formulation of the problem, Yang-Mills theory is supplemented with 
Higgs fields, and the response of a medium to the presence of non-Abelian 
charges ( at the macroscopic level) is described by two phenomenological 
constants in a phenomenological Lagrangian. An attempt of this kind has 
recently been made in , where the approximation of gluostatics developed 
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earlier for the case of trivial vacuum was used to describe nonrel- 

ativistic heavy particles carrying non-Abelian charges. It was shown that 
the equations of gluostatics with Higgs fields are consistent and that these 
equations describe Yang-Mills fields generated by the charges themselves and 
the response of the vacuum to the presence of these charges. Some general 
properties of solutions to these equations and the possible scenarios of charge 
confinement have also been discussed. 

This article reports on some results obtained from a numerical analysis of 
the equations of gluostatics involving Higgs fields. The results of numerical 
calculations are used to construct a simple approximation for the potential 
interaction between non-Abelian charges. By means of a direct substitution, 
it is shown that there are vortex-type excitations in the Yang- Mills- Higgs 
system and that these excitations are identical to Abrikosov vortices. Owing 
to this, it turns out that the numerical results for the tension coefficient that 
were obtained in p| for a funnel-shaped potential of interaction between mag- 
netic charges in a superconductor directly apply to chromomagnetic charges 
in a nontrivial vacuum. 



1 Chromoelectric charges 

In the model in question, the Lagrangian density has the form 

L = -\g^''G,, - rA, - \d^xD,X - \{X^ - F^}\ (1) 

where and x ci-re the triplets of, respectively, Yang-MIlls and Higgs fields 
( hereafter, the analysis is restricted to the case of the SU{2) group); is the 
density of the external-source currents; G^y = d^Ay — d^A^ + gA^ x Ay is the 
strength tensor of the gluon field ( the cross denotes vector product in isotopic 
space); and the covariant derivative is defined as D^ifs = d^ip + gA^ x ip. 
The equations for the Yang-Mills and Higgs fields are obviously given by 

D^G^,y + gDyX X X = Ju, 

(2) 

D^^D,x-Hx'-F'}x = 0. 

In the gluostatic approximation, we retain only the zeroth component of 
the 4-current. We have 

f = p = Pi6{x-xi) + p2 6{x - X2), 



2 



where P, the vector-charges and coordinates of the particles 

involved. Further, the vector-charges of the particles represent a convenient 
basis — Pi,P2, Pz = Pi^P2 — in which we expand solutions to equations (0). 
The zeroth component Aq = (p oi the gluon field then appears as a linear 
combination of the vector-charges of the particles, 

Cp{x, t) = ipi{x)P^{t) + ^2{x)P2{t), (3) 

while the vector field is proportional to the third component of the three 
basis vectors: 

Aix,t)=aix)P,{t). (4) 

The potentials and the vector fields are functions of the spatial coordinate 
and particle coordinates: ipi{x\xi, X2), a{x\xi,X2). 
By virtue of the condition 

d''j^ + gA''xj^ = 0, 

which ensures the consistency of equations (g), the basis generally proves to 
be rotating in isotopic space in the direction of the vector Q=^iPi+ V^2-P2, 
with frequency whereas the functions (pi=(pi{x2) , and 'P2=f2{xi) 
determine the values of the potentials and the points where the charges reside. 

Nontrivial solutions to the system of equations (|) were studied in [|| for 
the Higgs field represented as a linear combination of the vector-charges of 
the particles 

X{x,t) = xi{x)Pi{t) + X2{x)P2{t). (5) 

Following factorization of the contributions of the resulting functions a , 
ip^ = , V'sll , and = Wxi ^ X2 ||, the system of equations of 
gluostatics assumes the form 

DD^ + Jx}Cx = 5, 

DDx + g''{^Jx}C^ = KxJCx 

VxVxa- gj^ + gj^ = 0, 
j^ = <!>JD(!>, j^ = xJDx. 

where the column $ is the difference of the columns ip and ^, ^ = \\ ^1, ^2 



F'}x, 



(6) 
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Dki = V5fci + gaCki { k,l = 1,2) is the covariant derivative, 

6'^ = \\6{x — Xi), 6{x — X2)\\, and C and J are the 2x2 matrices 

1 

-1 ■ 

The parentheses in the above expression denote scalar products of vector- 
charges in isotopic space, and 6 is the delta function that describes the charge 
source of intensity equal to unit. 

A more general formulation assumes the use of additional components 
that have nonzero projections on other basis vectors as well. In contrast 
to electrodynamics, non-Abelian theory involves an extra degree of freedom 
associated with the choice of the relative orientation of the Yang-Mills and 
Higgs fields in isotopic space. The resulting system of equations is rather 
cumbersome, but the regularities in the changes that its solutions suffer as 
the result of this generalization may be traced by considering, instead of (^, 
the limiting case in which the Higgs field have a nonzero projection only onto 
the third vector of the basis: 

X{x,t)=x{x)p3{t). (7) 

Instead of the system of equations (H), we then have 

- gy^<i> = 6, 

Ax-g'<^\ = \{x'-F'}x, (8) 
VxVx a - gj^ = 0, = ^JD^. 

Let us supplement this system of equations with boundary conditions. 
The delta sources on the righthand sides of (^) and are eliminated by 
isolating the Coulomb term in the solutions for the field (p as 

Experience gained from the calculations and the analysis performed previ- 
ously for the equations of gluostatics in the trivial vacuum reveal that, at 
large values of the coupling constant g, we must approximate pointlike source 
more accurately, taking into account induced charge density [0. However, 
for the parameter values considered below, the above-type superposition in- 
volving the Coulomb solution and simple boundary conditions at the points 
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where the charges are located are quite sufficient for achieving a reasonably 
high accuracy of the calculations and a reasonably fast convergence of iter- 
ations. The quantities ^ were treated as parameters in one version of the 
calculations and were determined in solving numerically the problem with 
free boundary conditions of the second kind at the points where the charges 



are located 



0, 



dip2 



X=X2 



= 0. 

le rotating basis is that the boundary con- 
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An advantageous feature of t 
ditions at spatial infinity ( when the calculation is performed in a large box) 
are specified in a very simple form: |g — 0, ttp^z |g 0. In this case, the 
Higgs field must approach the vacuum value = F'^, whence it follows, in 
particular, that, if the charges are equal in magnitude (|Pi| = IP2I = -P), we 
have Xi |g = X2 |g = Xas = F/{2P cos(6'/2)), where 6 is the angle between 
the vector-charges in isotopic space ( this solution is obviously singular at 
6 = n). In the case specified by equation (J^), the boundary condition at 
infinity has the form x\g = F- The formulation of the problem admits the 
possibility of imposing additional ffist- and second-order boundary condi- 
tions on the field x that are analogous to the above boundary conditions for 
the fields 9? at the points where the charges reside. 

The potential of interaction between the two charges is determined by 
integrating over space the symmetrized field-energy density 



T = £, + 



mi 



= \e^JCE^ + ^lp,\^JxY + \e^JCE^ + -^{x' - F']\ (9) 

where E^ = E^ = Dx, and H = \/xa. 

Integration by parts makes it possible to recast the expression for the 
energy density Se of the gluoelectric field into the form 

= -yjC5 - €p^{^^Jx}{^Jx} - iP^{VJE^}a + 

2 (10) 



which is more convenient for numerical calculations. 

Let us ffist consider the case in which the Higgs field is specified by 
equation (^) and analyze solutions to the system of equations (||) . The energy 
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of the gluomagnetic field and the contributions to the energy that come from 
the five terms in (|^) are presented in the table as functions of the distance 
between the charges [ after integration over the space of energy densities with 
our numerical solution to the system of equations @]. For the parameters, 
we choose the following characteristic values: g = 1, P = 1, A = l, F = 1, 
and 9 = n — n/lO ( in connection with the choice of 9, it should be noted that, 
in numerical calculations, the vector charges cannot be taken to be strictly 
parallel because Xas is singular in this configuration). The calculations were 
performed for the free boundary conditions at the points where the charges 
are located. If the boundary conditions of the first order are used and if 
the parameter values are chosen as above, the minimum of energy can be 

achieved by varying ^] as a result, we arrive at the same solution. It should 
also be noted that the singular contributions that describe self-interaction 
were eliminated in the same way as in electrodynamics. 

From the table, we can see that the fourth and fifth terms represent the 
largest corrections to the first term, but they are of opposite signs and com- 
pensate each other. The second and third terms and the contribution of the 
gluomagnetic energy are negligibly small. Following a natural regularization 
of divergences, the first term is expressed in terms of ^ and takes the form 



Figure 1 displays the potential Vint as a function of the distance r between 
the charges for F = 0, 1, 2, and 3 ( A, g, P, and 9 are specified as above). For 
the trivial condensate F = 0, we naturally arrive at the Coulomb interaction 
represented in Fig. 1 by points. As F is increased, the curves are shifted 
below, this shift being equidistant at large r. 

In discussing the data presented in the table, we have already noted that, 
owing to cancellation of the leading corrections, the interaction potential is 
determined primarily only by the values of the fields cp at the points where 
the charges reside. Thus, the detailed distribution of the vector and Higgs 
fields may prove not very important, and even a rough approximation to the 
numerical solution will yield a reasonable result for the interaction potential. 
Indeed, let us assume that the vector field is negligibly small in the system 
of equations (1), and let us approximate the higgs field in the entire space 
by its asymptotic value. The equations for the fields ip then become trivial. 
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and the required solution is found straightforwardly ( the angle 9 drops from 
final expressions). As a result, we obtain 



^2 



Y ^-gF\x-xi\ \ _ Q-gF\x-x2\ 
SttIx — Xil SttIx — a;2| ' 

(12) 

_ ^-9F\x-xi\ \ Q~gF\x-X2\ 



■in \x — Xi\ ovr \x — X2I 



In the difference — ^21 the Coulomb components cancel out completely. 
If this were not the case, the fourth term in expression would guarantee a 
linear growth of the potential ( this provides a convenient tool for checking 
the accuracy of our calculations) ^ . 

The required constants 'y? are now given by 

^1=^2= p - 13 

QTX \Xi — X2I OTT 

The interaction potential can then be approximated as 

v,,,,.mi\l±JL!lL^A. ,14) 

Stt [ r J 

The following comment is in order. Only for A ^ 1 could we hope that 
the above crude approximation is reasonable. This is the limit in which the 
Higgs fields differ noticeable from their asymptotic values in the region of 
dimension r < A~^/^ and make a negligibly small contribution to the total 
energy integral. But owing to cancellation of the significant corrections, this 
approximation proves valid down to A ~ 1, which is a surprising result. 

Figure 2 shows the interaction potential as a function of the coupling 
constant (7 at a fixed distance between the particles ( r = 5). On the curve 
presented in this figure, we can find the d value corresponding to the onset 
of a nonlinear regime ( it is quite conceivable, however, that this restriction 
is peculiar to the computational algorithm that we used and has nothing to 
do with real physics). 

It was indicated in that the mode of solutions that is characterized 
by the suppression of the vector field \a\ =0, = must manifest itself 
with increasing A. Figure 3 shows the energy of the gluomagnetic field as a 
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function of A. We can see that the vector field in the system is suppressed as 
the superconducting properties of the medium are enhanced. 

It only remains for us to consider the case described by the system of 
equations which corresponds to the Higgs field proportional to the third 
basis vector. It can easy be seen that the charges are now coupled only 
by short-range forces because, in this field configuration, the gluoelectric 
fields are screened at large distances, the screening factor being e'^^"^. For 
the same reason, the generation of gluomagnetic fields is insignificant. As a 
result, the charges do not interact at large distances, so that, instead of 
we approximately have 

Vint = — • (15) 

Att r 

From the above results, it follows that solutions of the type (J^) are ener- 
getically favorable in the case of attraction between the particles [ (-P1-P2) < 
0], whereas the configuration (|^ is preferable in the case of repulsion. All 
these solutions are stable in classical field theory. Solutions are determined 
( "strengthted" ) by the boundary conditions, which specify the choice of the 
transverse component of the Higgs field lying in the plane spanned by the 
vector- charges of the particles. An effective potential that is expected to arise 
in quantum theory will contain the above solutions with a factor of continual 
integration. 

A simple comparison of the potential (0) with the potential used in the 
model of heavy quarkonia ( see f^), 

<y 

V{r) = + pr + Vo, 

r 

where a = 0.27, (3 = 0.25 GeV^ and Vq = -0.76 GeV ( r is measured in 
GeV~^), leads to the following estimates of the constants: 

g ^ 2.6, F ^ 1.1, ( — ^ 0.27, - — = 0.76). 

8 TT 8 TT 

At the same time, potential models yield a much lower value for the coupling 
constant g 

giJ/ij) = 1.55, giT) ^ 1.42, a,(J/^) ^ 0.19, a,(T) ^ 0.16. 

Conceivable, this is indicative of the importance of loop corrections. The 
relationship between the results presented in this study and those produced 
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by the potential models of heavy quarkonia will be studied elsewhere in 
greater detail. 

In this section, we have discussed the case in which the non-Abelian 
charges are placed in a strong gluonic superconductor. For the sake of com- 
pleteness , we present an approximation for the case in which the vacuum 
reduces to a trivial one ( A — >■ 0, F — > 0). For not overly large coupling 
constants ( g'^/{Ar[) < 1), the interaction potential that takes into account 
the contribution of the gluomagnetic field has the form 

{P,P,) + 7P| 6 - 7rV2 

Vint = , 7 = ^ 77 • 16 

Air r An Ibn 

2 Chromomagnetic charges 

It was noted in that the system of equations (J^), which corresponds to the 
approximation of gluostatics, can also be used to describe magnetic charges. 
We will illustrate this statement by means of an explicit substitution. In the 
absence of chromoelectric charges, we begin by eliminating singular delta- 
function sources from the system of equations and seek a particular solution 

of the form $ = ( in the case under study, the column is considered as a 
free parameter). 

The resulting system of equations, 

DDx = \{xJCx-F^}x, 

(17) 

VxVx a + gj^ = 0, = xJDx, 

describes the expulsion of a gluomagnetic field from a Higgs condensate just 
in the same way as this occurs in the Abelian theory of superconductivity. 
From the symmetry properties of the system of equations and from the fact 
that the eigenvalues of the matrix C are complex conjugate to each other 
( /^i,2 = T^l-Psl), it follows that a Higgs doublet can be described by one 
complexvalued function. The equations of the theory are then equivalent to 
the Ginzburg-Landau system of equations. Explicitly, this can be demon- 
strated by means of the substitutions 
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X2 



1 - e'^ ~ TTe^ 



2|A 



Instead of two real-valued functions, transformations (|T8]) define a complex- 
valued function ip ( here, ip is the complex conjugate of ip)- In addition, 

TP 

we introduce the scale transformations V —>■ qF V and a —>■ a, which 

reduce the equations for the function ip and the vector field a to the canonical 
form 

(19) 

VxVxa - + ia)^p - {V - ia)ip} = 0, 

where k = \. 

9 

It is well known [|lO] that the system of equations (0) has solutions with 
a quantized magnetic flux. In the initial variables, these solutions correspond 
to Abrikosov vortices with the magnetic flux 

f adl = Qn = n = ±1, ±2, . . . , 

where / is a closed contour that circumvents the vortex un the unperturbed 
condensate. Formally, the flux of the gluomagnetic field is then given by 

i Adl ^ 



J 9 IP3I 

The gluoelectric field identically vanishes 

E = PiD^i + P2D<^2 = , 

and so does the gauge-invariant 't Hooft tensor of the field strength; that is, 
we have 

F,. = TGI, - ^ r (^mX X D^xT = 0, 
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where x = The energy density is concentrated near the vortex core 

and is given by 

E = g'F'S, £: = l(Vxa)2 + l|(V-m)^P + |{|^p-ir. (20) 

At K = 1/ \/2, the integral of energy is quantized and has the hnear density 
J Ends = 71 F'^n. In the Abehan case, solutions with large \n\ are stable for 
K < 0. In the case being considered, the gluomagnetic field does not 

necessarily disturb a gluonic superconductor, because there is an additional 
degree of freedom associated with the relative orientation of the Yang-Mills 
and Higgs fields in isotopic space. Thus, solutions of the Abrikosov-vortex 
type are metastable, but they may lie comparatively far from the instabil- 
ity region. Solutions of this type have already been used in the theory of 
electroweak interaction and are referred to as W and Z strings |]12| . 

Unfortunately, magnetic charges separated by large distances can be de- 
scribed only by solving numerically equations ([T9|) with the input singular 
Dirac potential 

^ g j z-d z + d \ 

° A'K p\[p^+{z- rf)2]V2 [p2 ^ ^ ^)2]l/2 J ' y ) 

where e^p is a unit vector, p and z are the radial coordinates, and rkd are 
the points at the z axis where the magnetic poles reside. The results of such 
numerical investigations are presented in p[ . among other things, the tension 
coefficient for a funnel-like potential was calculated in these studies as a 



function of the parameters of the Ginzburg-Landau potential. Formulas (18), 
which establish relationship between the Abelian and non-Abelian models, 
can be used to extend these results to the case of chromomagnetic charges 
in a gluonic superconductor. 



Conclusion 

The problem of two nonrelativistic chromoelectric and chromomagnetic char- 
ges in a Higgs vacuum with the properties of a gluonic superconductor has 
been considered within classical field theory. It has been shown that there is a 
direct analogy with the Abelian theory of superconductivity, where magnetic 
charges are confined, while electric charges are not ( by charge confinement. 
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we mean here a linear growth of the energy of particle interaction with dis- 
tance). It has been found that Abehan modes do not exhaust all possible 
states of chromoelectric charges; for these, there exist states whose energy 
is less than the Coulomb ( Yukawa) energy. It is conceivable that the in- 
clusion of loop corrections may reconcile the results presented in this study 
with available data on the potential of heavy quarkonia, thereby ensuring 
confinement of charges. 

A simple representation is given above for a certain class of solutions to 
the Yang-Mills-Higgs equations that describes non-Abelian Abrikosov vor- 
tices. The solutions obtained here can be used in constructing models of 
elementary particles, which can then be considered either as closed vortex 
hues in the bulk of a Higgs condensate or as droplets of a dual condensate 
with a quantized gluoelectric field frozen into them. In the latter case, the 
source and outflux of the gluoelectric field from the droplet surface can be 
interpreted quark- ant iquark pair. 

Considering the dual pattern, we arrive, for the dual fields ( potentials) 
and charges, at the scenario in which the fates of chromoelectric and chromo- 
magnetic charges are interchanged: the former are confined, while the latter 
are free. 
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Gluoelectric and gluomagnetic contributions to the energy 
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